The seed magnetic field generated during recombination 
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Nonlinear dynamics creates vortical currents when the tight-couphng approximation between 
photons and baryons breaks down around the time of recombination. This generates a magnetic field 
at second order in cosmological perturbations, whose power spectrum is fixed by standard physics, 
without the need for any ad hoc assumptions. We present the fully general relativistic calculation 
of the magnetic power spectrum, including the effects of metric perturbations, second-order velocity 
and photon anisotropic stress, thus generalizing and correcting previous results. We also show that 
significant magnetogenesis continues to occur after recombination. The power spectrum ^/WPb 
decays as fc^ on large scales, and grows as fe°'® on small scales, down to the limit of our numerical 
computations, ~ IMpc. On cluster scales, the created field has strength ~ 3 x 10~^^ Gauss. 

PACS numbers: 98.80.-k,98.58.Ay 



I. INTRODUCTION 

Evidence is growing for magnetic fields on larger and 
larger scales in the Universe (see e.g. the reviews [Hill). 
In galaxies, the fields have strength of order ^Gauss, 
ordered on scales ~ 1 — lOkpc. Fields of strength 
~ 1 — on scales ^ 0.1 — IMpc have been de- 

tected in galaxy clusters, and there is evidence of mag- 
netic fields in superclusters. Recently, new evidence has 
been presented for intergalactic magnetic fields: high en- 
ergy gamma-rays from distant sources can initiate elec- 
tromagnetic pair cascades when interacting with the ex- 
tragalactic photon background; the charged component 
of the cascades will be deflected by magnetic fields, affect- 
ing the images of the sources. Using observations from 
FERMI, a lower bound of order 10~^® G has been claimed 
for the strength of fields in the filaments and voids of the 
cosmic web 

The origin of these fields is still unclear (see e.g. [^-Q)- 
They could have been generated via astrophysical pro- 
cesses during the nonlinear collapse stage of structure 
formation. There remain unresolved difficulties in ex- 
plaining how these astrophysical seed fields lead to fields 
of the observed strength and coherence scales. Alterna- 
tively, the fields could be primordial seed fields - created 
in the very early Universe, during inflation, or during 
subsequent phase transitions. In principle inflation can 
generate fields on all scales - but unknown physics must 
be invoked to achieve non-minimal coupling of the elec- 
tromagnetic field. The electroweak and QCD transitions 
can only produce fields on very small scales, up to the 
Hubble radius at magnetogenesis (and their amplitude is 
strongly constrained by their gravitational wave produc- 
tion before nucleosynthesis jlG|). 

Primordial magnetogenesis also takes place in the cos- 
mic plasma after particle/anti-particle annihilation. This 
avoids the problem of exotic physics that faces inflation- 
ary magnetogenesis - standard Maxwell theory and stan- 
dard cosmological perturbations in the cosmic plasma in- 



evitably lead to magnetic fields. It also avoids the small 
coherence scale problem facing electroweak and QCD 
fields. However, the problem is the weakness of the fields, 
since this effect occurs at second and higher order in cos- 
mological perturbations. 

The key question is how weak is the field and how 
does it vary with scale? Differing qualitative estimates 
of the field strength have been given by [TT| - [l6j . The 
power spectrum was first numerically computed by 1171 . 
which differs significantly from ours. More recently, [l8| 
presented a power spectrum that is closer to our result. 
We discuss below the differences between previous results 
and ours. Our analysis is the first complete general rela- 
tivistic computation of the power spectrum, taking into 
account all effects. 

Our result is shown in Fig. [1] The power spectrum 
behaves as 



fc"^ fc < fc, 
fc > fc, 
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On cluster scales the comoving field strength is 

SiMpc-Sx 10-29 G. (2) 



Thus the field generated around recombination is too 
weak to act as a seed for the observed field strength of or- 
der ^G. Adiabatic contraction of the magnetic flux lines 
during nonlinear collapse of structures provides an en- 
hancement of ~ 10"^, while the nonlinear dynamo mech- 
anism has an amplification factor ^ 10^ (with many re- 
maining uncertainties). Note that hydrodynamical and 
turbulence effects during nonlinear collapse themselves 
generate a field of order 10"^° G - which is also too small 
to account for the observed galactic and cluster fields Q . 

The field ^ is also too weak to imprint detectable ef- 
fects on the CMB. Nevertheless it is a real property of the 
standard cosmological model, and may have some impact 
on early structure formation during the 'dark ages' if it 
is the only primordial field. (See e.g. [l^ [2^ for the role 
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FIG. 1: Left: Magnetic field spectrum today. Right: Comoving magnetic field strength today at a given scale. 



of magnetic fields in structure formation during the dark 
ages.) 

As shown below, the magnetic field is given by 



3e 
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(3) 
(4) 



where $, are first-order metric perturbations, Au^..^ — 

vl — v!y is the photon-baryon velocity difference, and 0] is 
the photon quadrupole moment, from anisotropic stress. 
This leads to three types of source terms for magnetoge- 
nesis: 



$W_vI/(l)}A^;(;'] 



(5) 



The first source term is second-order, while the other two 
are quadratic in first-order quantities. The contributions 
of the source terms to the power spectrum are shown in 
Fig. [5](left). 

Our paper builds on the physical ana,lysis of nonlinear 
plasma dynamics presented in [isl - ITsl . [2l| - [23| . The key 
features of the dynamics are as follows. 

• The electric field ensures that the proton-electron 
relative velocity is always strongly suppressed in 
comparison with the photon-electron relative ve- 
locity - even at high energies when the Compton 
interaction is stronger than the Coulomb interac- 
tion. 

• Vorticity induced in the electron fiuid is thus trans- 
ferred almost entirely to the protons, and the 
baryon vorticity evolution is determined by the 
two-fluid dynamics of photons and baryons, which 
is very close to the equations of CMB dynamics. 
We use the second-order Boltzmann code of (23| . 

• The limit V(. — Vj-^0 and Wp — Wo — ^ is not equiv- 
alent to setting Vp = Vc = Vj in the momentum 



exchange equations, and the limit must be taken 
consistently. 

• At first order, cosmological vector perturbations 
are zero after infiation, in the standard model. 
Magnctogenesis requires vortical currents, and 
these can therefore only be generated at second or- 
der, via mode-mode coupling of first-order scalar 
perturbations. This remains true even in the pres- 
ence of topological defects, which are active sources 
for vector perturbations: at first order, the vector 
perturbations induced by the defects cannot break 
vorticity conservation in the cosmic plasma (2^ . 



• On large scales there is some cancellation amongst 
the source terms in ([S]) (this is evident from Fig. 
[5]). Neglecting any of the effects can thus lead to 
unreliable results. 



• The magnetic field continues to be created after re- 
combination, due to the residual nonzero ionization 
fraction. If the numerical integration is stopped at 
recombination, then the comoving field is under- 
estimated by a factor ~ 10 (see Fig. [5]). 



The plan of the paper is as follows. In the next section 
we review and clarify the magnetic and electric field gen- 
eration beyond the tight-coupling limit. In Sec. IIIIl we 
detail the numerical integration of the differential evolu- 
tion equations at second order in cosmological perturba- 
tions that we perform in order to solve for the magnetic 
field spectrum. We also provide analytical insight into 
the time and scale behaviors of the numerical results. 
We compare our results with previous work in Sec. IIVI 
Details of some calculations are given in the Appendices. 
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II. UNDERSTANDING THE ORIGIN OF THE 
MAGNETIC FIELD 

A. Interactions in the cosmic plasma 

The stress-energy tensor of a species s satisfies 



(6) 



The momentum transfer rates are given by 



, Aw^ EE 

3C ' pC 
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C^^ ^ qsns{E'' + e^,.<B") , s = e, p . 



(12) 
(13) 



/3-^,(14) 
mp 



(15) 



where C!f^.{— —C^^) encodes all the effects of interactions 
with species r. Relative to observers with 4- velocity u^, 
the energy density transfer rate is — it^C^, and the mo- 
mentum density transfer rate is = h'^C^^^ where the 
projector is h^^ = 5"^ + u^u^ . 

The Euler equation for a species s is given in general 

by 



(7) 



The radiation energy density p^, the quadrupole of the 
radiation temperature anisotropy 6^;^, and the number 
densities rig, are as measured by observers. In the 
rest frame w^, the electrons and protons are well approx- 
imated by pressure- free matter, T^^ = pl'^^*^u^u^, where 
Ps^*^' is the rest- frame density measured by u^. In the 
frame, there is effective pressure, momentum density and 
anisotropic stress: T/"' = psuf^u" + Psh'"' + 2qi'"u''^ +77^" , 
where I2lt. 



The kinematics of are described by decomposing its 
covariant derivative as [2ll ^ 



(8) 



where is the volume expansion, tr^i, is the projected (i.e. 
orthogonal to u^^), symmetric and tracefree shear, cj^,y is 
the projected antisymmetric vorticity, and ii^^ = u'^V^Uf^ 
is the projected acceleration. The vorticity vector is de- 
fined as 



AW 



(9) 



where the totally antisymmetric tensor is defined by 
£0123 = \f~~9- (Note that our sign convention for w^i/ 
and definition of recover the Newtonian limit, and 
differ from [2l|,[2|.) 

In the period of interest, from the end of particle/anti- 
particle annihilation up to now (T-,, < 500keV, z < 
2 X 10^), the relevant species are protons, electrons, pho- 
tons, and when recombination occurs, hydrogen atoms. 
Neutrinos affect only the background dynamics and the 
gravitational potentials in the Einstein equations. The 
Faraday tensor of the electromagnetic field defines elec- 
tric and magnetic fields measured by observers: 



1 



(10) 



Protons and electrons couple to the electromagnetic field 
through the term C^p = F^^jsi where s = p, e and is 
the electric 4-current. Then Vt.T^'' = - Y^^F^jl- We 
have = QsUsU^, where qs is the particle charge, Ug is 
the number density (in the rest frame) and the 4- velocity 
of species s is 



[l-v!)-'^"- (11) 



Here 7s is the relative velocity of s measured by u^. 
Maxwell's equations are given in Appendix Rl 



2 rest 
7s Ps 
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Ps = ^v!p, , (16) 
Ps(«-^t'.'/»^''). (17) 



The Thomson cross-section is ctt = Sna'^ / {3ml), and the 
Coulomb interaction is governed by the electrical resis- 
tivity 



Vc 



Tre^^TTio In A 
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^ . (18) 



where A is the Coulomb logarithm. On cosmological time 
scale s the m agnetic field diffuses below a length scale 
^ yVc/Ms ^ 100 AU, so that diffusion can safely be 
ignored 18|. The characteristic time scales for electrons 
interacting via the Coulomb and Thomson interactions 
are 



rrir 



rc 



TT 



20 sec / 1 + z 



-3/2 



0-TP7 



5 X 10«sec 
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l + z 

103 
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,(19) 



(20) 



where ric is the number density of free electrons and Xc 
is the fraction of free electrons. We used n^o + 'Iho — 
3 X 10~^cm~3 [2^. The time scale which characterizes 
the evolution of the plasma can be taken as 



,(^) 



min{rs(z), TiMpc{z)} 
1 



mm 



^/H{z)aTnc{z) (1 + z)ki Mp, 



TzAl 1-(21) 

1 + z)«:iMpc J 



Here rg is the Silk damping time and 1 Mpc is taken as 
the minimum comoving scale on which we can trust a 
second-order pcrturbative analysis up to redshift z = 0. 
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Electric field 



The Eulcr equation 
velocities is given by [21 



7]) for the proton and electron 



msUs V. 



(22) 



where s, r = p, e and 



u^i^s H Kg D^Tis + D^Wg^ ) vl: 



(23) 



ng 4 

Tig 3 



The covariant spatial derivative is defined in (jA6p . 
The first term on the right of ((231) describes not only the 
evolution due to the expansion of the universe which con- 
serves the particles, but also the evolution of the number 
density due to recombination which does not conserve 
the particles when hydrogen atoms are formed around 
recombination. 

Prom now on we expand in perturbations around a 
Friedmann background, up to second order. The metric 
in Poisson gauge is 

ds^ = [-(1 + 2$)dr/2 + 2Sidx'dri + (1 - 21')da;2](24) 

where Si is a vector perturbation (9*S'i = 0) and en- 
ters only at second order. Perturbed quantities arc ex- 
panded according to X = X + X'-^'> + X^'^'> . Only the 
first order of scalar perturbations $ and ^ will enter 
the evolution equation of the magnetic field, so we omit 
the superscripts for them. The explicit form of the term 
i)^^ + + in (HU is then given by ((CT|) . with 

WS = — Cg. 

We set ric ~ rip = n, since we find that the final ex- 
pression of the resulting electric field is not affected by 
n-c — np, in agreement with [2^ . 

In order to obtain a dynamical equation for the velocity 
difference Av^^ 



— v^, we use ([22]) to obtain 



men 



+ (5{C^^^+C^^) . (25) 



The Lorentz force term in (|15p has been neglected since 
it is higher order. We define the baryon velocity as the 
velocity of the centre of mass of the charged particles; 
then 



(mp + mc)w,^ = TOpw(^ -t- rUcV^ , 



(26) 



1 + /3 



vt; = iT - ^-tA<, . (27) 



'"^ 1 + 13 



In principle, the baryon velocity can be different from 
the velocity of hydrogen, i.e. of electrons and protons 
rccombined, but thermal collision ensure that hydrogen 
atoms follow closely the electrons and protons. 



Using (I25|)-(l27l) and the explicit forms (IT2l)-(IT4]) of the 
collision terms, we obtain 



Av^^ + AK^,) = (I + P)eE'' - (1 + P)e^n7jcAv^, 
(l-/33)(A„^^ + ^ex) 



l + /3^ 
1 + P 



[av^, + le^^Av;, 



(28) 



We show below that the B^jAup^ term can be neglected, 
since it is higher order. 

Equation ([28)) shows that an electric field can be gener- 
ated by nonzero velocity differences Avpe and At;^b- The 
Maxwell equation (|A2[) shows that then can be gener- 
ated, provided that is transverse. We will show that 
the generated electric field keeps electrons and protons 
more bound together and therefore leads to a decrease in 
Aiipo, which becomes negligible compared to Av^c- 
Neglecting third order terms, the Maxwell equation 
[3]) can be rewritten in terms of the velocity difference 



Aw^o as 



A<^ = —( curl B" - E"^ --eE''+ a^''' eJ] , (29) 
en V 3 / 



where we used (|A5|) . 

In order to estimate the magnitudes of the various con- 
tributions in the stationary regime, we expand all evolv- 
ing quantities in frequency space: 



dcj7\/''(a;,a;)e''^'' , 



(30) 



where the mode has characteristic oscillation fre- 
quency uj ~ Tp~o- I'^ terms of the characteristic timescales 
pQ]) and dSO]), we find from ^ and ^ that 



(i+/^)+of^+4^+i 



V T, 



Am, 



3 Tcvo 

VcTc 



(l + /3) + ofi^^) 



TcvoTT 

curlB^ 



3eTT 



-(1-/3^) 



?7C,cff7'cvo' 

5 



(31) 



where we used A/Vpc = 0{Avpc), and we defined ^ 

4(1 + /34) rc 



VC,eS = VC 



1 



3(l+/3)2rT 



(32) 



Given the hierarchy of the different timescales involved 
in (|3T|) . it follows that the largest contribution to the 
resulting electric field is given by the velocity difference 



Ai 



b7 ' 



This can be seen in Pig. [21 where we plot the 



different ratios of typical timescales that enter in pip . 
Specifically, all the plotted ratios are always well below 
unity for the period of interest, from very large redshift 
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FIG. 2: Evolution with redshift of different ratios be- 
tween characteristic times that arise in (|3ip . compared with 
unity (thick black line): r/cTc/ir^vo) (thin solid), rjc/rcvo 
(thin dashed), 77ct"c/(tovoTt) (dotted), 77cTc/(??c,GfiT-cvo) (dot- 
dashed) and rjc,eS /tbvo (thick dashed). (The jumps in the 
curves occur at reionization.) 



until today, even accounting for recombination around 
z ~ 1080. This allows us to vifrite 



ilc,cs curls'" 

4mc 1-/3^ 
~ 3eTT 1 + /3 



(33) 



In order to compute the final magnetic field produced by 
such an electric field, we consider the curl of the electric 
field, governed by Maxwell's equation (|A2[) . In frequency 
space 



-rjccS curl curl 
^ 3ct^ 1 + /3 



curl AD^^ + • (34) 



Remembering that the magnetic field is divergence free, 
we can compare the first two terms of the above equation. 
Their ratio in Fourier space is of order (Tovo'7c,off^^)^^ — 
''cvo/?7c,cff • Therefore, on all scales of interest, we can 
conclude that the contribution of the rye. oft curl iJ^ term 
in ((33|) is negligible compared to the last term. 

The above considerations remain valid once we ap- 
proach recombination time, as long as the residual frac- 
tion of free electrons Xc is not too small. This is to ensure 
that the approximations of the ratios of time scales made 
to obtain (|33|) remain valid. This is indeed the case, and 
it can be checked from Fig . ^ since x^, ^ 1Q~^ — 10~^ 
after last scattering (26l . |27| until reionization. 

We are therefore left with the following expression for 
the electric field produced by the tiny velocity difference 
between electrons and protons: 



1 - Ap^CJT 

1 + 13 3e 



(35) 



It is important to note that this expression does not con- 
tain the number density of free electrons Ue- Therefore 
the electric field produced by this mechanism before re- 
combination is still present after last scattering (see also 
j26| ) and can in principle continue to generate a magnetic 
field after recombination. 

We can now also finally prove that 



Aw'' < Aw'' 



(36) 



Using ((35|) and (|34|) without the rjceS term, and in the 
Maxwell equation (|29p leads to an estimation of the order 
of magnitude of velocity differences: 



TcvoTT 



PI 



Aw'' 



(37) 



The order of magnitude of the ratio Av!^^ 



/Av^^ is shown 



in Fig. [5] and remains well below unity for all relevant 
times, even when Coulomb scattering becomes less effi- 
cient than Compton scattering, that is for z > 10^. 
It also follows from (|37| and (fT3)) that we can rewrite 
as 



(38) 



where we neglect terms of order /? and where here the 
baryon index b encompasses protons, electrons and hy- 
drogen atoms. 

As a conclusion of this section, we stress again that 
when we assume that electrons and photons are tightly 
coupled, as was originally considered in [28j . then the 
electrons and protons are even more tightly coupled by 
the electromagnetic field which is generated, so that the 
electrons and protons can still be considered, from the 
point of view of photons, as a single fluid of baryons. 
As a consequence, taking Av^^ —J- at early times has 
to be performed consistently by keeping Av^^ <^ Aw^^ 
when taking the limit. For the tight-coupled limit, this 
is crucial, since it corresponds exactly to the limit Ve = 

= Up = 0, and the collision terms cannot be evaluated 
directly from their expressions P^ - P^ . 



C. Local inertial frame (tetrad) 

It is convenient to express all quantities in a local in- 
ertial frame, deflned by an orthonormal tetrad Ba (a = 
0,1,2,3): 



ah 



(39) 



The tetrad indices are distinguished from general coor- 
dinate indices by underlining, and «, j, fc • • • = 1, 2, 3. We 
choose a comoving tetrad, so that eg is the fundamen- 
tal observer 4- velocity: eg'' = u^- In the background, 
eo'' ~ ~ (a~^,0). The perturbed tetrad is given in 
Appendix |B] Derivatives along the tetrad vectors are 
defined by 



(40) 
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Covariant derivatives in the tetrad frame are computed 
using the afBne connections given in Appendix [B] 

Tetrads make the physical meaning of all nonscalar 
quantities more transparent. In linear perturbation the- 
ory, it is common practice to decompose perturbed quan- 
tities in a background tetrad. For instance the veloc- 
ity is often decomposed as u^^^^ = a"^?;^^^^, together with 

ul^^ = civ^^\ which means implicitly that v['^^ = JyW^j^j. 



where here, as in ([38|) . the baryon index b encompasses 
electrons, protons and hydrogen atoms. Finally, note 
that the value of the magnetic field depends of course 
on the observer. Its value in the baryon frame is related 
to its value (|42p in the fundamental frame by 



i£k 771 



(47) 



Thus v'f^^^ coincides with v^^-^ = e-^Uj-^^. Introducing 
tetrads is the natural generalization of this standard pro- 
cedure when considering higher order perturbations, and 
this has already been used for example to decompose ve- 
locities [i^, llOl • The nonlinear evolution of the distribu- 
tion of photons is well suited to computation in a tetrad 
frame 1231. 



E. Numerical computation 



In order to solve the evolution equation for the mag- 
netic field, we need to solve the Boltzmann hierarchy for 
baryons and photons, to compute the source of the elec- 
tric field in (|35|) . The basic idea is to decompose the 
directional dependence of radiation in the local inertial 
frame into multipoles: 



D. Magnetic field 



The Maxwell equation ()A2|) becomes in the tetrad basis 

doia^B^) = -ah^e [(1 + * - Ek] , (41) 

Equivalently we can use derivatives in the coordinate ba- 
sis: 



d^fc 

(27r)3/2 

47r Ni/2 



^er(fe)G,™(fc,a;,n)(48) 



GeUk,x,n) = e''=--V^™(n^) . (49) 

We suppress the time dependence for convenience. 

Terms quadratic in first order perturbations appear as 
convolutions, and we introduce the notation 



.^B^y . -a\m, [(1 + <!> - *) E,] , (42) ,c{f,mk) . J SUk. + - k) Mk,) Mk,) , 



where we have used the fact that the electric field is at 
least a first order quantity, and the magnetic field a sec- 
ond order quantity. The gravitational potentials in this 
expression occur only at first order. Equation (|42|) is 
compatible with [? ], which can be seen via Ek = eu'Ei. 
To obtain (Hi}, we need 



(50) 

A Fourier mode qi is decomposed on the helicity basis of 
the background spacctime as 



= = + g(-)e(_) + 9(0)6(0) , (51) 

= H^lil) ■ (52) 



{cuA EY- = ei^Ve_Ek = e^^de [(1 - ^)Ek\ , (43) 

which uses the affine connections up to first order given 
in Appendix [Bl Also, 



e^^e^-'^u.E^^e^ueEk 



which follows from 



e^Epdk^ 



(44) 



{u^W^u,)eY' = {WoeK)e^'' = = d,^ . (45) 



In addition, we omitted terms like ^e—diEk_ and 
'^e—de_Ek in deriving (|42|) . since the electric field con- 
tributes only at first order - and at this order, it is curl- 
free. For the same reason, we can also replace by 
a~^dt. 

In summary, magnetogenesis is governed by (j42[) and 
i.e. 



The background helicity basis vectors f!.(h)i with helicity 
/i = 0, ± are defined in [23|. The azimuthal direction 
h = corresponds to scalar perturbations and is aligned 
with the total Fourier mode, i.e. e(o) = fc, while h = 
± correspond to vector perturbations. At first order, 
when the mode is aligned with the azimuthal direction 
since q — k, there arc only scalar perturbations. For 
vector quantities like the electric field, we need to use a 
helicity basis e(;j) on the perturbed spacetime, and this 
is built by the identification of e(/i) with e{h), i-e. e^^^ = 
^\h)- Vector quantities like the electric field Ei arc then 
expanded as 

X- = ^(+)4+) + + X 



(0)6(0) 



(53) 
(54) 



a'BA = 



e(7ic -I- nii)x, 



(1 + $ - *) Cj 



b7 



In this basis, the Maxwell equation becomes (ex- 
plicitly giving the perturbative order of quantities) 



- —f ; ^ — ^iikO 

e(rio + n]i)Xc — 



(l + $-*)V,T"^ ,(46) 



(fc) + /C{ [$(1) - } (fc)] . (55) 
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Wc projected (|42)) along ef^''* and used 



(±) 



±ke 



(±) 



(±) 



(56) 



Note that there are only contributions from h — ± and 
we thus recover that scalar perturbations cannot generate 
a magnetic field and vortical perturbations are required 
to source the magnetic field. Using the multipole decom- 
position of (j35p . and neglecting /? ^ f , we obtain finally, 



(fe) 



= ±ka' 



3e 



{[5W+<l>W-vl/(i)]F(2)}(fc) 



±fca 



2 4crTP7 



3e 



5^±)(fc) + 5f'(fc) + 4±)(fc)l (57) 



where 



and (J-y — Spj/p-y. Also, 
k(/i,0) = ^/{4~Wj, K{h,±l) = - 



(58) 



{2±h){3±h) 



(59) 



The last equality in ([57)) defines the contribution of each 
line above: 5*^^^ is the purely second order contribution 
from l/*^^'; is the S^V contribution and Sg^'' is the 
02fb contribution. 

Although ^^^^(fc) vanishes at first order since there 

are no vector perturbations, V^'_^|(A;i) and V^'_^|(fe2) do 
not vanish in general, since the modes ki and k2 are not 
necessarily aligned with the azimuthal direction k = k/k. 
We first need to obtain their expression when the modes 
fei or k2 are aligned with the azimuthal direction, and 
then we perform a rotation of the azimuthal direction 

In order to explicitly take into account the symmetry 
of the convolution products in (j57p . we can symmetrize 
the source terms. At first order there are only scalar 
perturbations, and all first order tensorial quantities are 



gradients of scalar functions, so that X, 



(1) 



X 



(1) 



di-^ . . . di^^X^-^\ Most of the source terms are of the form 
e—de (XYk) = e—de (XdkY), and once projected along 

fib)* 

e] they contribute to the generation of the magnetic 
field proportionally to 



^ / J^W^'^^^^ t^^^ ~ ^^^^'^^ " X{q)Y{k - g)].(60) 

Here X and Y denote S^.V^^^Vb,^,"^. 

This symmetrization, which is always possible, shows 
that for these types of terms, the configurations of 



{k,ki,k2) with ki = k2 will not contribute in the convo- 
lution. Only couplings from a quadrupolar quantity to 
gradient terms, which are of the type 

e^9, (Xk'djY) ^ S^dt {dkd^Xd.Y) , (61) 



as in the last line of ([57| . can have contributions to 
the convolution coming from configurations with fci = 
k2. The generated magnetic field is thus severely sup- 
pressed at early times for these configurations since 
the quadrupole of radiation is suppressed in the tight- 
coupling regime. 



III. NUMERICAL RESULTS 
A. Transfer functions 

In order to obtain the final magnetic field spectrum 
produced via this mechanism, wc integrate numerically 
the evolution equations for cosmological perturbations up 
to second order, since we have to take into account even 
the behavior of the second order velocity difference be- 
twecn baryons and photons V^f^-j {k , rj) . We use through- 
out the cosmological parameters of WMAP7 [3l| . 

For a variable X, the first order transfer function is 
X(i)(fc,77) = A'(i)(fc,7?)$in(fc), where is the gravita- 
tional potential deep in the radiation era. Because of 
statistical isotropy, the first order transfer function de- 
pends only on the magnitude of the Fourier mode and 
not on its direction. This is however only strictly true 
for multipoles like 9™ and V(^h) defined from non-scalar 

quantities if the azimuthal direction is aligned with fc, 
and considering only scalar perturbations at first order 
the contributions for h ^ vanish. However, when us- 
ing these first order transfer functions in the quadratic 
terms of the second order equations, we must rotate these 
multipoles according to the angles between fei, k2 and k. 
This is to ensure that the multipoles remain defined with 
respect to the total momentum k [23|. 

The second order transfer function X^^^{ki,k2,'ri) is 
defined by 

A(2)(fc, Tj)=IC {X^^\k,,k2, 77)$i„(fel)$in(fc2)} (k) . 

(62) 

Without loss of generality we enforce A'^^^ (fci, ^2, 77) = 
X^'^'f {k2, ki,ri) in numerical calculations. The transfer 
functions of the first and second order quantities needed 
in the source terms are obtained by a joint solution of 
the Boltzmann equation (for photons and neutrinos), the 
conservation and Euler equations (for baryons and cold 
dark matter) and the Einstein equations (for metric per- 
turbations). They are found numerically using the same 
techniques as in |32| . 

The transfer function of the magnetic field can be split 
into the different contributions of the S^^^ sources de- 
fined in (|57p . The transfer functions of these contribu- 
tions are related to the transfer functions of the sources 
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through 

(63) 

and this is how we obtain the complete time behavior of 
the magnetic field. A crucial point that will turn out to 
have important consequences is that the final rcdshift for 
numerical integration should be taken after the recom- 
bination epoch. The electric field that results from the 
small electron-proton velocity difference and that gives 
rise to a magnetic field is still present after last scatter- 
ing, when the fraction of free electrons Xc is tiny but still 
does not completely vanish (see also [26j). 

In order to compute the equal time correlation func- 
tions of the magnetic field, we need the power spectrum 
of the initial potential, defined by 



{^UkmM))^Kk-q)P{k). 



(64) 



If the source terms are Gaussian random variables, we 
can apply Wick's theorem, and the contributions of the 
two polarizations ft. = ± add up quadratically; 

{B{k,ri)B*{k\r,)) 

-^^^^)3^/d^9P(<z)P(|fc-g|)x 

1 (9, fe - q, r/) p + {q,k-q, ??)S(*+) {k - q, q, r/) | 

= ^^Jy'^ J d'qP{q)P{\k - g|)|g(+)(g,fc - q,v)\' 
= 5l{k-k')PB{k,r)), (65) 

where S(-i-) = ^i^f±)- In the last line we have defined 
the power spectrum of the magnetic field Pb- Its value 
today is plotted in Fig. [5l 

In order to have a deeper analytical understanding of 
the resulting magnetic field spectrum, we study each con- 
tribution Si independently. There are cross correlations 
in (|65p . but our aim is to assess the relative importance 
of the different contributions; the Pj' are defined by re- 



placing with Sj^j in (|65 



B. J^-Aub^ contribution 

The velocity difference between baryons and photons is 
severely suppressed in the tight-coupling limit relative to 
other perturbations like 5^] we expand this tiny velocity 
difference in terms of the expansion parameter k/r' <^ 1, 
where r' = ricaTa is the derivative of the optical depth 
for Thomson scattering. At first order in k/r', in the 
radiation-dominated background on supcr-Hubblc scales. 



vlo)\k,v) 



R- 

T 

const . 



k J vL 



(66) 
(67) 



Using R = 3pb/(4p'y) oc a, 1/r' oc a ^ and a oc 77, we get 

4+'(|fc-q|,g,77)cxg(+)(<7'-|fc-q|3)^ . (68) 

Then ([55)1 gives the early-time and large-scale behaviour 
of IS^^y and the resulting magnetic field power spectrum 
behaves as 

P|^(fc,7y) cx e j <:\^q\q(+)\^P{q)P{\k-q\) 

'7oq 

For a scale-invariant initial power spectrum, P{q) oc g^"^, 

Pi^Xk,7j)^X'PlHk,ri), (70) 

as can be seen just by a change of variable in the integral 
of dini). In [13 it is found that Pi^{Xk,r]) = X^P^^k,T]). 
The disagreement appears to arise since |18| infer the de- 
pendence on k from the q ^ k contribution to the integral 
in (|69p - but the main contribution to that integral are 
also limited to g < fc given the argument at the end of 
section Hi El We finally find that for the S2 source term, 
the power spectrum of the magnetic field behaves as 



k^P^^k,r]) cx k^^ 



(71) 



This behaviour in k and 77 at early times when the mode is 
still super-Hubble, is confirmed by numerical integration, 
as is evident from Fig. |3] (left). 



C. 02«b contribution 

Similar analytical arguments apply to the magnetic 
field generated by the source 53. The tight coupling ex- 
pansion of the source is 

02(fcj'7) « -j'^'o « — ' «b(o)(fc,'y) « fc??, (72) 

T-' Vcq 

in a radiation background on supcr-Hubble scales. This 
implies that the S3 contribution to the magnetic field 
power spectrum behaves as 



k^P^''{k,ri) cx k"^ — 

Vcq 



(73) 



It has the same k dependence as ([7T|) but a different ?/ 
dependence. The analytical form is verified by the nu- 
merical output shown in Fig. [3] (right). 



D. Aujj^ contribution 

For the purely second order part 5*1, the only way to 
assess its contribution is to consider the tight coupling 
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5^.Av contribution for k/yfceq=.5, 1,2,5, 10,20 ©2. v,, contribution for k/A:eq=.5, 1,2,5, 10,20 




Logiol^/°eq] Log,o[a/ne,] 

FIG. 3: Left: Magnetic field spectrum Pg^{k, 77) from only the 5*2 contribution in (|57p . for different k/kcq, with values increasing 
from bottom to top. Right: Magnetic field spectrum Pg^{k,rj) from only the S3 contribution in (|57p . 



expansion of the evolution equation for the vorticity of 
baryons. Indeed, we need to evaluate first the total con- 
tribution Si at lowest order in tight-coupling, and the 
detail of this derivation is given in Appendix IC 21 It 
follows that behaves as {k/T'){k'ri)'^ oc k^ff/rjl^, 

which implies that for the total magnetic field 



^k^PBik.-q) cx fc*-^. (74) 

'7cq 

This behaviour is confirmed by numerical integration, 
as shown in Fig. 2] (right). Since 52 oc k^ / Tj'^q, 
S3 oc k^rj^/ricq, and ^^^V^ /Vcq^ obtain that 

iSi oc k^ri'^/rjaq- Thus ^3 contributes to the magnetic 
field power spectrum as 

\A^^FiM)«fc'— , (75) 
which is verified in Fig. 2] (left). 

E. Magnetic power spectrum 

From these plots it is evident that the magnetic field 
is still generated after recombination. This is the reason 
that it is important, to set the final time of integration 
after recombination, since the largest contribution comes 
from this last period of generation. Indeed, before reach- 
ing the usual 'final' stage where the magnetic field is no 
longer sourced but only redshifts with time {B oc a~^), 
we observe a bump in the resulting magnetic field spec- 
trum, corresponding to the recombination time. This 
should be interpreted as an increase in magnetic field 
generation due to decoupling of photons and baryons. 

In the decoupling regime the fluid of photons and 
baryons is no longer equivalent to a perfect fluid. The de- 
parture from tight coupling may be interpreted via non- 
adiabatic pressure perturbations, which can source the 



total vorticity [TBI, Issl - issj . It is not a priori evident that 
this could lead to an increase in the magnetic field gen- 
eration. On the one hand, the total vorticity is sourced 
when interactions between baryons and photons are less 
efficient, but on the other hand, there is less vorticity ex- 
change between photons and baryons since the collisions 
are less efficient. In the ideal limit where the decou- 
pling is complete, the vorticity of photons and baryons 
is adiabatically evolving according to (jC6p . whereas the 
total vorticity is sourced by the gradients in the total 
non-adiabatic pressure. This is possible because the vor- 
ticities of the different fluids do not add up linearly to 
give the total vorticity as can been seen from (jC4p . 

However, when decoupling occurs, we observe that 
there is in fact an increased generation of magnetic field 
in that phase, and this essentially comes from the factor 
Xo in (|46p . i.e. from the fact that the magnetic field is gen- 
erated via the residual ionized fraction. More precisely, 
the generation of the magnetic field is proportional to 
9[jV^Tj^j,j/a;o and not only to d[jV ^T^ so even when 
V^r^j, around decoupling, "S/ ^T^ j./ Xc can still have 
sizeable values. This last significant stage of magnetic 
field generation is counterbalanced and finally stopped 
by the redshifting of photon energy density {p^ ex. a~*). 
It can be seen from ([57]) that the background radiation 
energy density controls the efficiency of the total mag- 
netic field production after recombination. 

The power spectrum of the magnetic field is shown in 
Fig. [5] (left). The behaviour on large scales (oc fc^) is ex- 
plained above. The behaviour on small scales is complex, 
since it depends mainly on the generation between hori- 
zon crossing time and Silk damping time. During that 
period, the analysis which we restricted to super-Hubble 
scales does not apply - and the adiabatic redshifting does 
not apply either, since the magnetic field continues to be 
generated. For k ^ kgq, a reasonable linear approxima- 
tion is log (^/k^Ps) oc 0.5 log k. 
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Av'^^ contribution for k/iteq=.5,l,2,5,10,20 



All contributions for k/A:en=.5, 1,2,5, 10,20 



10-^3 



10- 




10- 



100 



o 

s 

10-=' 




Log,ola/fleq] 



0.1 0.5 1.0 5.0 10.0 

Logiol^/feq] 



50.0 100.0 



FIG. 4: Left: Magnetic field spectrum Pj^(fc, ry) from only the 5*1 contribution in (|57p . for different k/k^q, with values increasing 
from bottom to top. Right: Magnetic field spectrum PB{k,rf) for all contributions. 
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FIG. 5: Left: Magnetic field spectrum today (solid). Contributions from the different sources in ((5]) are distinguished: second 
order velocity term Si (dot-dashed), quadratic term S2 in velocity and density (dashed), quadratic term Sz in anisotropic stress 
and velocity (dotted). Right: Comoving magnetic field strength at a given scale at times 1 + 2 = 1, 10, 100, 1000 corresponding 
respectively to solid, dashed, dotted and dot-dashed lines. (Dashed and solid lines cannot be distinguished). 



F. Magnetic amplitude 

The magnetic field amplitude smoothed over a comov- 
ing scale A is 

= dfcfc2pB(fc)cxp(^-^j, (76) 

where the normalization volume is = 
/ d3yexp[-yV(2A2)] = X^{2T:f^^. Note that the 
integral is insensitive to the upper cutoff, which may be 
taken to infinity, since A ^ Adamp- The magnetic field 
strength is shown in Fig. [5] (right). 

The field strength at 10 Mpc is approximately 10^^^ 



Gauss and three times as much on cluster scales 1 Mpc. 
Given the slope of the spectrum, this is expected to grow 
to larger values for smaller scales. Our numerical inte- 
gration does not allow us to investigate smaller scales 
since the numerical integration time increases dramati- 
cally with fcniax- In addition, the results become unre- 
liable on small scales where density perturbations have 
become nonlinear by z = 0. On the comoving scale of 
the Hubble radius at equality, the strength is ^ 10""^° G. 



G. Frame dependence 

At early times when photons and baryons arc tightly 
coupled, the magnetic field measured in the baryon- 
photon fluid is not generated at lowest order in the tight 
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coupling expansion. This is shown in Appendix ID 21 
Only higher orders in the tight-coupling expansion con- 
tribute to magnetogenesis. However, since most of the 
magnetic field production occurs when the tight-coupling 
expansion breaks down around recombination, this sup- 
pression is only relevant at early times, before recombina- 
tion, and for modes which remain for the longest time in 
the tight-coupled regime, i.e. for large scales. Therefore 
the difference between the magnetic field in the funda- 
mental frame and in the baryon frame decreases, and 
they are nearly equal today, as shown in Fig. [S] This 
shows that at 1 -f z = 1000 there is a suppression for 
large scales in the baryon frame, but today there is no 
more suppression since most of the magnetic field has 
been generated around recombination time. 



{l+zfBx(G) at 1-Hz=l,1000 




5 10 50 100 500 

A in Mpc 

FIG. 6: Magnetic field strength at a given scale as measured 
in the fundamental frame at 1 -I- z = 1 (continuous) and 1 -I- 
z = 1000 (dashed), and as measured in the baryon frame at 
1 + 2 = 1 (dotted) and 1 + z = 1000 (dot-dashed). Dotted 
and continuous lines cannot be distinguished. 



IV. DISCUSSION AND COMPARISON WITH 
PREVIOUS RESULTS 

Our approach is the first complete analysis of magne- 
togenesis around recombination, in the sense that it does 
not neglect any term in the second order equation for 
the generation of the magnetic field - previous work has 
omitted at least one of the terms. Therefore our results 
will necessarily differ from existing partial results and we 
discuss briefly how some of these differences arise. 

Two general points can be highlighted: 

• Numerical computation is essential to obtain the 
magnetic power spectrum - and even for a reliable 
estimate of the magnetic field strength. For ex- 
ample, [ll,[ll use similar analytical methods and 
incorporate the same source terms, but the two esti- 
mated field strengths on the recombination Hubble 



scale differ by orders of magnitude. A full numer- 
ical integration is needed, especially to take into 
account all orders in the tight-coupling expansion. 
This was initiated by [l^, and we have built on 
their work. 

• Neglecting any of the source terms for magnetogen- 
esis not only leads to inaccurate predictions - it also 
misses the fact the separate source terms do not 
simply add up linearly. The total of the different 
contributions is suppressed in the tight-coupling 
regime on large scales by a factor {kij)^: the details 
are given in Appendix ID II As a consequence, dis- 
carding some terms implies that this suppression in 
the tight-coupling regime is neglected - which leads 
to an overestimate of the magnetic field generated. 
This is especially critical for the largest scales where 
tight-coupling is valid at the latest times. 

In [ll|, [l^ [l3| the anisotropic stress contribution, 5,3 in 
(|57p. and the second-order velocity contribution. Si, are 
neglected. It is apparent from the power spectrum plot 
in Fig. [5] that both of these contributions are substantial 
and cannot be neglected for a reliable prediction of the 
magnetic field. In addition, these references omit the 
scalar metric perturbations. Metric perturbations and 
the second order velocity are included in [TB - flTj , but the 
anisotropic stress is neglected. 

In [3| the anisotropic stress is included, but the second 
order velocity contribution is neglected. In addition to 
this difference from our work, we find a different time and 
momentum dependence for the large-scale and early-time 
behaviour of the ^2 and S3 contributions. We then find 
\/WPb fx while they find oc fc*"/^. 

The first numerical prediction of the magnetic power 
spectrum was given by [l3| , neglecting anisotropic stress 
but including second order velocity. However, our power 
spectrum is significantly different from theirs. Part of 
the difference is due to anisotropic stress, but there is 
a further difference arising from the treatment of veloci- 
ties. The evolution equation for the magnetic field can be 
given by (|46p . It is true that in the tight-coupled regime 
(see Appendix ICl for details), the velocities of electrons, 
protons and photons can be approximated to be equal. 
However, it is erroneous to use e— cJfV^T^^, = to esti- 
mate the vorticity evolution. Indeed, in order to cancel 
the collision term when taking the tight-coupling limit, 
we have to consider a combination which uses the action 
reaction law and for which the collision terms do not ap- 
pear. It is given by the total fiuid vorticity conservation 
equation: 

e^c),V^T;, + e^deV^T^, = . (77) 

In the tight-coupled limit, the fiuid of baryons and the 
fluid of radiation exchange vorticity, essentially because 
the dilution of their energy density is different, and this 
exchange of vorticity is then required to maintain equal 
velocities at all times. In [13] it is implicitly assumed that 
Cti^ can be neglected because the velocity of electrons 
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is assumed to be close to that of photons. However, as 
we discussed in Sec. |lTl the hmit has to be consistent 
with (|37l) . and this colhsion term is precisely responsible 
for the vorticity exchange between photons and electrons, 
and thus between photons and baryons - and it cannot 
be ignored. The vorticity evolution in the tight-coupling 
limit should be computed using (jD6[) . i.e. by substitut- 
ing the tight-coupling solution of velocities and energy 
densities perturbations in (|46)) . 

In [l^ it is shown that there can be no generation of 
magnetic field in the photon frame at strictly less than 
the first order in tight coupling (if there is no initial 
vorticity). Note that what we call first order in tight- 
coupling (see also [sB] ) is called second order in tight cou- 
pling by [l^ [3| • In our case, we focus on C/^^ , whereas 
they focus on {k/T')C!^^ where r' is the interaction rate 
and k/r' is the parameter of the tight-coupling expan- 
sion. The result of [Tsj is compatible with our results in 
Appendix ID 21 since in the tight-coupled regime the pho- 
ton frame is the baryon frame. Thus the magnetic field 
in the photon frame will be generated only starting from 
the next order, i.e. at first order in the tight-coupling 
expansion. Our numerical approach does not rely on a 
tight-coupling expansion since we integrate the full sys- 
tem of equations, and in that sense we consider neces- 
sarily the full tight-coupling expansion in our computa- 
tion. We checked numerically that at early times, when 
photon-baryon coupling is efficient, the magnetic field in 
the baryon frame is severely suppressed compared to the 
magnetic field in the fundamental frame. 

V. CONCLUSION 

We have performed for the first time a full numeri- 
cal computation of the seed magnetic field generated by 
nonlinear dynamics, taking into account all general rel- 
ativistic effects and all source terms. We discussed the 
range of applicability of the mechanism on cosmological 
scales and concluded that the generation of the magnetic 
field is directly related to the Compton drag by photons 
on baryons. Even in the tight coupling regime, photons 
exchange vorticity with baryons and the magnetic field is 
created. Since the electric field that sources the magnetic 
field does not depend on the fraction of free electrons, 
the magnetic field is still generated after recombination, 
given that there is a relic fraction of charged particles, 
and we find that the largest production takes place in 
this final stage. 

Our results are summarized in Fig. [1] The power 
spectrum (left plot) behaves as 

On cluster scales the comoving field strength is (right 
plot) 

Si Mpc~3x 10-29 G. (79) 
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Appendix A: McLxwell's equations 

Maxwell's equations V[AF^i.] = and V^.F'"' = in 
a general spacetime take the form [l^ [25j | 

D^B^ = -c^^F^, ^^E^' =u:,,B^ + Q, (Al) 

= -curl E^-e^.^ii-'E^ (A2) 
+ loE^ - - uj^,)E-' 
= curl Bf, + e^,Au"S^ - , (A3) 

where E^^,B^^ are defined by ([10]). Here the total 4- 
current is = j^' + and it is split as 

J'^ = ^?^i^ + J^ g=-u^f, J''=Kf, (A4) 

where g, are the charge density and current measured 
by u'^ observers. By pT|) . 

g = e(7pnp - 7ono), = e(7pnpW^ - 7cnoW^). (A5) 

The derivative is the projected covariant derivative 
and it defines a covariant curl [2l|, [2^ : 

D^./ = /i;:V,/, D^^"^ = hlK^^S\ (A6) 
curls'^ = e'^'^^D^S'A. (A7) 

We work in Gaussian units so that the fine structure 
constant is a / {An) = 1/137.036 and the magnetic 
field strength is measured in Gauss. 
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Appendix B: Tetrads 

The tetrad basis is given up to second order in scalar 
perturbations by 
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(Bl) 
(B2) 
(B3) 
(B4) 



This choice of tetrad is discussed in [23[ (see also [29|, |37|, 
113). The covariant derivative of a tensor in the tetrad 
basis is given by 

VaX,^ = e^a^X/ - + , (B5) 

where indices are lowered and raised with 'i]ab and 77—. 
The afhne connections in the background are 



and the perturbed forms are 



"ifeO — (Jik , rt — — 

— a ~ a 



"OOi ~ "oiO ^ ' "Oifc " " ' 



(1) 



(1) 



(B6) 

(B7) 
(B8) 
(B9) 



Appendix C: Euler and vorticity equations 

1. Euler equation 

For a perfect fluid with equation of state Wg = Palps 
and speed of sound Cg = dPs/dps, the term on the left 
of the Euler equation ^ is given to second order in the 
tetrad basis by [2^ . i39i l : 



Ps{l + Ws) 
l + c2 

+-; ^ 

1 +Ws 



uf + (1 - 3c^,)HK + 



1 + w, 

,2\„,s 



OA + 



- m'u\ 

2\ 



<' + (1 - -ic'Mul + — 

- 1 + Wi 



1 + 



- 



3H(1 + m;s)2 



(CI) 



2. Vorticity evolution 

The vorticity tensor of species s is 



(C2) 



and the vorticity vector is given by ([9]). In the tetrad 
basis, up to second order, 



(C3) 
(C4) 



The evolution of the vorticity is deduced from ([7|) and 
()Cip . For a non-interacting perfect fluid, up to second 
order [H, HI 



Ps{l+Ws) 



,\' + (2-3Cs2)Ha;,^fc=0. (C5) 



This can be recast as 

[ps(l + Ws)a^w|]' = 0. 
For an interacting fluid, 
(2 - 3c2)Hw|, 



1 + Ws 



(C6) 



}.(C7) 



Appendix D: Magnetogenesis in tight-coupling 

1. Magnetic field in fundamental frame 

In the case where there are only interactions between 
baryons and photons, C^^ + C^^^ = 0, and 



(Dl) 



In the tight-coupled limit where the interaction rate be- 
comes very high, photons and baryons behave like a single 
fluid, with 



1 



1 



Iff = 



3 + 4i? ' ''^ 3(1 + i?) ' '~ Ap^ 
The energy density contrasts at first order are 



i? ^ ^ . (D2) 



5(1) ^ 



(D3) 



The velocities of baryons and photons arc the same in 
this regime 



b 7 f Z 7 2 



By (leg) and (jPTj) . 



0~wlfc +H(2-3c^)w*fc 



[pf(l + Wf)aVfci' 
Pf (1 -I- Wi)a'^ 



(D4) 



(D5) 
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This can be used to infer the source term for magne- 
togenesis in pS)) . In the tight-coupled regime, 9[jV^T^^ 
can be estimated by using (|D4[) and (jD3p in the baryon 
version of (|C1[) . Then, subtracting d^i^ ^T^^-^ — 0, we 
obtain 



-^%M/9fc]<5f|. (D6) 



From pS)) it then follows that in the tight-coupled 
regime, the evolution of the magnetic field is given by 



2. Magnetic field in baryon frame 



From (|T7|) we obtain 



Bl-B!- 



^iik h TP 



aeXo (1 4- W{) 



VgdkSf , 



(D9) 



where the second equality holds in the tight-coupled 
regime. Using the first order version of the Euler equation 
(|C1[) for the plasma, i.e. with V^Tj^- = 0, and using also 
the first order evolution equation for the plasma density 
contrast. 
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,2 

■f Mk 



m 

l+Wi 



(1 



1 



1 + Wl 



d[i<^dk]5f \ ,(D7) 



V 1 + wf / 



3*' - a-w^ 



we deduce that in the tight-coupled regime 



a'Bt 



{"'My 



(DIO) 
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where we used ph = (wp + mc){nc + nn) — inp{nc + nn). 
Note that Sc'^H = dln[pb/(/Ob + 4/3p^)]/d77. Since the 
vorticity in the tight-coupled plasma obeys (|D5p , the first 
term on the right hand side of (jD7p . which is linear in 
the vorticity, can only source the magnetic field if there 
is initially vorticity in the plasma. This is the term re- 
sponsible for the Harrison mechanism [23, . All other 
terms which are quadratic can source the magnetic field 
even if there is no initial vorticity. 

However, on large scales in the radiation era there is 
a suppression of the total contribution of these quadratic 
terms. From the large-scale radiation era relations at first 
order. 



2'Hdi4 '- 



-2$. 



(D8) 



it follows that at lowest order the quadratic terms are 
estimated by diXdjY ~ di^dj^. Hence the quadratic 
source terms are suppressed by a factor (kr])'^, since 
at lowest order all contributions are of the type ^ 
d[i^dj]^ = 0. This implies that y/k^PB{k,ri) cx 
kWhlld^ that is Y,,^S^ cx k^ij^T]^^. 



At early times in the radiation era we have a;e — 1, and 
then we obtain a conservation equation up to second or- 
der: 



-Bt 



(D12) 



This is precisely the Harrison mechanism, but up to sec- 
ond order. 

In the tight-coupled regime, in the plasma frame, the 
magnetic field can only be generated if there is initial vor- 
ticity, i.e. through the Harrison mechanism. We recover 
here the results in [l^ [l^ . The magnetic field measured 
in a different frame is only due to the contribution of 
the electric field to this change of frame. In the fun- 
damental frame, this contribution in the tight-coupled 
regime is given by the second and third lines of (|D7p . 
Note that the electric field is generated at first order in 
cosmological perturbations even in the lowest order of 
the tight-coupling approximation and even in the plasma 
frame. 
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